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Abstract 

We present a matrix technique to obtain the spectrum and the analytical index of some 
elliptic operators defined on compact Riemannian manifolds. The method uses matrix 
representations of the derivative which yield exact values for the derivative of a trigono- 
metric polynomial. These matrices can be used to find the exact spectrum of an elliptic 
operator in particular cases and in general, to give insight into the properties of the so- 
lution of the spectral problem. As examples, the analytical index and the eigenvalues of 
the Dirac operator on the torus and on the sphere are obtained and as an application of 
this technique, the spectrum of the Dirac-Kahler operator on the sphere is explored. 
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1 Introduction 



The problem of finding the spectrum and eigenstates of elliptic operators on compact 
Riemannian manifolds is of considerable interest in physical theories, from quantum field 
theory to gravitation and cosmology. The spectral problem of the Dirac operator on the 
sphere, for instance, is related to confinement of quarks; it has been studied since many 
years ago and has been solved on the n-dimensional sphere [H [21 El H]. The spectral 
problem of the so-called Dirac-Kahler operator acting on differential forms, postulated by 
Ivanenko and Landau in 1928 as an operator to describe properties of fermionic fields and 
connected with the Dirac equation by Kahler in 19620 is also physically relevant since it 
has been used to describe differential forms as spinors at the time that it exhibits bosonic- 
like solutions [H El [7J [8] , giving an alternative framework for lattice theory [9] [TUl [Til [T2] . 
In a sequel of papers (p3]-[T6], and references therein), a nonlocal formulation of lattice 
theory has been introduced. In some simple cases this technique yields exact results, i.e., 
the numerical output can be interpolated to obtain the functions that solve exactly the 
problem. In this method, the discretization of a differential operator is obtained by the 
substitution of the derivatives appearing in the operator by matrix representations of the 
derivative which arise naturally in the context of interpolation of functions. Once the 
finite-dimensional operator is constructed, it can be solved by standard techniques. 
The aim of this paper is to show the potential use of this lattice formulation by obtaining 
the spectrum and analytical index of the Dirac equation on the torus and on the sphere 
and exploring the eigenvalue problem of the Dirac-Kahler operator on these surfaces. 
We solve the Dirac operator on the sphere and find that besides the known fermionic 
spectral solution with eigenvalues ±/, I = 1, 2, . . ., there are present a bosonic-like solution 
of the spectral problem with eigenvalues 0, ±3/2, ±5/2, . . . 

Concerning the Dirac-Kahler spectral problem, we find that the spectrum splits into two 
sets: one corresponding to the bosonic-like eigenvalues ±a//(/ + 1), I = 0, 1, . . ., and 
the other to a fermionic-like counterpart consisting of numbers approaching the integer 
eigenvalues of the Dirac operator on the sphere with double multiplicity, confirming that 
in general, the Dirac-Kahler operator cannot be decomposed in several (two in this case) 
independent Dirac equations on a curved manifold. 

It is worth to notice that the use of this technique allows to find analytically the Dirac 
spectrum on a multidimensional torus. 

2 Differentiation matrices 

In this section we present the main results on the discretization method based on the use 
of differentiation matrices. Proofs and further applications can be found in [151 HE]- 
Let us consider iV = 2n + 1 arbitrary points — 7r < x\ < x 2 < ■ ■ ■ < x 2n +i < 7r and the 

fact, the Dirac-Kahler equation is equivalent to several independent simultaneous Dirac equations 
in a flat manifold [5]. 
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N x N matrix defined by 
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where D and T are the matrices whose entries are given by 
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If f(x) is a trigonometric polynomial of degree at most n and denote the vector whose 
elements are the derivatives f^ k '( 
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1 N, then 



fM=D k f, k = 0,l, 



(3) 



If an even number of nodes, say iV = 2n, is selected instead of an odd number, this 
equation still holds but for antiperiodic trigonometric polynomials, i.e., functions of the 
form f(x) = exp(ix/2)g(x) where g(x) is a trigonometric polynomial of degree at most 
n — 1. 

If the periodic function f(x) G L 2 (— 7r,7r) is not a polynomial, a residual vector must 
be added to the right-hand side of 0- The norm of this vector approaches zero as the 
number of nodes is increased since f(x) can be expanded in a Fourier series. 
The differentiation matrix D takes the simple form 



0. 



J 



k. 
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2 sin §(.?-£;) 

when the nodes are chosen as the evenly spaced points 
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(2j-JV-l), j = l,2, 
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The generalization to the case of q variables straightforward. The space 

of functions where the differentiation matrices become the partial derivatives is the tensor 
product of the subspaces of trigonometric polynomials of degree at most in x k . Let 
iVfc be the number of nodes in the direction k. Then, the total number of nodes is 
N = Ylk=i and they are indexed by 



q-l 



jt + (j 2 - l)N x + (j 3 - l)mN 2 + ■ ■ ■ + (j q - 1) H N k , 



(6) 



k=i 
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where jj. = 1,2,..., N k , k = 1, 2, . . . , q. The N x N differentiation matrix representing 
d k /dx k , is now 

D fc = l Nq ® • ■ ■ ® l Wfc+1 ® D fe (8) 1^ ® • • • ® livx, (7) 
where k = 1,2, . . . ,q and it should be multiplied by the vector f whose entries are 

in order to obtain 9f(x^ l ,x^ 2 , ■ ■ ■ , x^)/dx k . 

The differentiation matrices ([7j) can be diagonalized simultaneously by the unitary and 
symmetric matrix 

F = F q (g) ■■■®F 1 . (8) 
Each Nk x iVfc matrix Ff, is a discrete Fourier transform with the property 

i,F),F) foFfc) j m = id^J 5j m . 

The components of Fk are given by 

(F k ) jm =^-e-^ d -^, (9) 

and the eigenvalues of D k are 

dg) = -(JV fc + l)/2 + m, (10) 

where j,m = 1, . . . , Nfr. 

3 Lattice Dirac operators 

In this section, we show how the eigenvalues, the eigenspinors and the analytical index 
of the free Dirac operator on the flat torus and on the sphere can be obtained from a 
discretized version of the Dirac operator constructed with differentiation matrices. 

3.1 The flat torus 

A two-dimensional flat torus is a Riemannian quotient manifold T\ = R 2 /A, where A is 
a discrete group generated by the translations on the directions Ai = (1,0), A2 = (a, b), 
where b > 0. There are four spin structures on Ta, each of them is characterized by 
(^1,1^2), where ^2 = 0,1 and has associated a spin bundle whose sections are spinor 
fields. A spinor field is a function / : IR 2 — > C 2 which is periodic or antiperiodic in the 
Afc-direction if \i k is or 1 respectively. 

In order to use differentiation matrices it is required a 27r-periodicity in the orthogonal 
coordinates x and y instead of the natural A-periodicity of the spinor fields. Since the 
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linear function ip(x, y) = (x + ay, by)/2ir in M. 2 yields a diffeomorphism between the torus 
M 2 /(27rZ) 2 and T\, we can use this change of variable to write the Dirac operator on T\ 
in terms of the coordinates x and y. Thus, the Dirac operator takes the form 

° D ~) (ID 
where 

In general, the spectrum of the Dirac operator depends on the choice of the spin structure. 
For T\ the spectrum of D is already known [T7]. If the spin structure of T\ is given by 
{(i>i, fa) the eigenvalues of (TTXj) have the form 



± 



> + f + + f)A£ (1.3) 



where A^ = (1, —a/b), \^ = (0, 1/6) and || • || stands for the euclidean norm in M. 2 . Each 
pair of integers m and n yields one eigenvalue of multiplicity K, where K is the number 
of possible combinations of m and n giving the same eigenvalue. 

The use of differentiation matrices yields a discrete problem with the same result as we 
show in the following. Since the eigenspinors are 27r-periodic functions in the fc-direction 
if /Xfc = and 27r-antiperiodic functions if = 1, we can use the differentiation matrix 
(Jlj) to yield the discretized form of the Dirac operator (TTT1) 



°-(£T). (") 

where 

D ± = (^±^)D lT iD 2 , (15) 

Dfc is given by ([7j). The spin structure will determine whether the number of nodes is 
odd or even. 

The eigenvalue problem of the discretized Dirac operator D can be solved by noting that 

det(D - Al) = det (Vl - (1 + ^)D 2 - ^D 2 , + ^D^^ , (16) 

where we have used the formula of Schur and the fact that commutes with Dj. This 
matrix can be diagonalized by the matrix F, where F is the two-dimensional discrete 
Fourier transform given by (jSJ) . Taking into account the structure of (|T6l) , we get that the 
eigenvalues of the discretized Dirac operator are of the form 
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where is one of the eigenvalues of given by fflQj) and || • || 2 stands for the usual 
hermitian norm in C 2 . If the spin structure is not twisted (the trivial structure), the 
eigenspinors are 27r-periodic trigonometric polynomials and must be an odd number. 
Therefore, S k > runs over Z when — > oo for k = l,...,n [cf. ( flOj) ]. If the spin 
structure is twisted in the fc-direction, the corresponding component of the eigenspinors 
is 27r-antiperiodic, must be now an even number and S k > has the form / + 1/2, where 
I is an integer. Thus, in the asymptotic limit N k — > oo, ( TT7|) coincides with . Note 
the relation between the spin structure (/xi,/x 2 ) and the parity of the set (Ni,N 2 ). The 
zero eigenvalue appears only in the trivial structure. 
The analytical index rip of the Dirac operator is defined by 

no = dimkerZ} 4 " — dimker.D~. 

Therefore, the finite-dimensional version of no is the difference 

n* D = dimker D + — dimker D~, (18) 

where D ± are the block matrices ffl5|) . Since the lattice Dirac operator is invertible, as 
the continuum operator, the index n* D is zero. However, we can compute explicitly n* D by 
counting the null eigenvalues of D ± . Again, these matrices can be diagonalized by F. In 
this way we find that the eigenvalues of D + have the form 

{l + f> ^ ~ \ d(2) (19) 

and those of D~ have the form 

(i_") d (i) + l d (2). (20) 
b b 

Therefore, both D + and D~ have null eigenvalues only in the case of the trivial structure, 
because only in that case, S 1,2 ^ can be zero. From (jTHJ) and f[2"Uj) . we can see that the 
number of null eigenvalues of D + is the same number of null eigenvalues of D~. Therefore 
n *D = n D, i-e., the analytical index of the discrete free Dirac operator on T\ is zero and 
this only occurs for the trivial spin structure. 

The above procedure can be straightforwardly generalized to higher-dimensional tori. 
3.2 The Dirac operator on the sphere 

This is another example of application of differentiation matrices yielding exact results. 
Let us consider the standard parametrization <fi : (0, it) x (0, 2n) — ► M. 3 of the unitary 
sphere S 2 given by 

4>{6,ip) = (sin 9 cos <p, sin 9 sin ip, cos 9) . 
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The hermitian Dirac operator takes the form 



D 



' 2— + -cotfl + csc#— \ 

\4 + r ot0 - csc ^ / 



(21) 



This problem has been solved in higher dimensions [HEJH]. In this case the eigenvalues are 
±1, ±2, . . . , ±n . . . with multiplicity 2n. The eigenspinors are given in terms of products of 
27r-antiperiodic trigonometric polynomials in 6 and (p. To derive a discrete Dirac operator 
we have to take into account that < 9 < n and therefore, the differentiation matrix 
Dq can not be constructed according to pj. However, we can use the general form (CQ) 
instead of PJ and choose Ng different points in (0,7r), or even better, take the unique 
solution < Q\ < 9 2 < ■ ■ ■ < 0n 6 < 7r of 

Ne ft —ft 1 

2oot2-^ = -iootfl i (22) 

Mi 

as nodes [16J. This selection yields a hermitian discrete Dirac operator D of the form (TH|) 
where now D ± are given by 

D± = il v ®D e ±D v ® esc 6, (23) 

where G is a diagonal matrix whose nonzero elements are the numbers 9k- The matrix 
corresponding to the term (cot^)/2 of (pTj) does not appear here because of (|22|) and 
therefore, Dg has null diagonal elements and it is a differentiation matrix for the operator 

1 d 



V sin 8 



acting on trigonometric polynomials. The differentiation matrix may be constructed 
according to (j4j) since there is no restriction on ip. Note that the number of nodes to 
be used in the construction of the differentiation matrices must be even integers in both 
directions. The eigenproblem Di/> n = \ n ij) n can be solved by the use of standard routines 
to produce the right results: if min(A^, N v ) is denoted by N, the exact eigenvalues 

±1,±2,-.. ,±k,--- ,±{N-2)/2 

appear among the eigenvalues of D with the right multiplicity 2k. In the limit N — > oo, 
the complete solution is recovered. The eigenvectors can be interpolated to give the correct 
eigenspinors but it is necessary to take into account that this is a degenerate eigenvalue 
problem. For example, let us take the eigenvalue Ai = 1. The corresponding eigenspinors 



e^/ 2 cos 0/2 \ . f-ie'^ 2 sin 9/2 

-xe^ 2 sin 6/2 J ' W> W - I e ^/2 cog Q /2 
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can be obtained as the linear combination of the eigensolutions of = ip\ . For 

Ne = N v = 10 we have that ip\(6,ip) can be obtained from axip^ + a2ip < { 2 \ where 
ai = 9.3217 — 3.1676i and = 0.8051 — 0.0499z. If ipi denotes the vector whose elements 
are the values ipi(0j, <p k ) ordered according to (jSJ), the error is given by 

||Vi - aiV! 1} " a3^{ 2) ||2 = 8.2 x 10" 14 . 

Concerning the index, the same simple argument used to show that the continuum Dirac 
operator has index zero works here: since the lattice Dirac operator is invertible, it has 
index zero. 

We end this section by noticing that there exist bosonic-like eigensolutions of the Dirac 
operator on the sphere. The components of this kind of solutions are 27r-periodic functions 
in both directions and they can be found through this technique by taking an odd number 
of nodes N$ and N v . The corresponding spectrum consists of zero, with multiplicity 2, 
and half-integer numbers ±3/2, ±5/2, . . . , ±(2n + l)/2 . . . with multiplicity In. 



4 The Dirac-Kahler operator on the sphere 



Let us consider the same parametrization of the unitary sphere given before and let us 
look for the solution of the Dirac-Kahler equation 



D DK (j) = rruj>, 



(24) 



where -Ddk — d — 8 and (f> is an inhomogeneous differential form on the open chart (0, n) x 
(0, 27r). Here, d is the exterior derivative and 5 the adjoint or interior derivative 5 = *d*, 
where * is the Hodge star operator. In the orthonormal basis (1, d9, sin 8dip, sin Odd A d(p) 
of inhomogeneous differential forms, the Dirac-Kahler operator takes the form 
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DK 
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de 
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csc z 9— 
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-oo +cotd 
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(25) 



/ 



To obtain a discrete Dirac-Kahler operator we proceed as before, with the same definitions 
of G, Dq and D v given in the last section. To find fermionic-like solutions of f !24p . i.e., 



27r-antiperiodic eigenfunctions in 6 and cp, we use an even number of nodes Ng and N v in 
the construction of the differentiation matrices. The numerical solution of the eigenvalue 
problem of the discrete Dirac-Kahler equation produces a two-fold set of real eigenvalues 
clustered around certain limit values, say m^, k = ±1, ±2, ■ • • , to which they converge as 
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Ng and N v grows. The number of eigenvalues clustered around rhk are twice the number 
corresponding to the multiplicity of eigenvalues of the Dirac equation, and the limit values 
rhk approach the eigenvalues of the Dirac operator on the sphere as k — > oo. The first six 
positive limit values rhk and their distance to the Dirac eigenvalues k, are shown in Table 1. 



Table 1 : The first six positive stabilized eigenvalues of the Dirac-Kahler 
operator on the sphere and their distance to the corresponding eigenval- 
ues of the Dirac operator on S 2 . 



k 




k — rrik 


1 


0.8660254 


0.1339746 


2 


1.9364916 


0.0635084 


3 


2.9580398 


0.0419602 


4 


3.9686269 


0.0313731 


5 


4.9749371 


0.0250629 


6 


5.9791303 


0.0208697 



It is worth to notice that using the basis \rshi8(l±i sin 8 d8 A dip), Vcsc 8(d8±i sin 8dip), 
the Dirac-Kahler operator takes the anti-diagonal form 

where D is the Dirac operator on the sphere given in (12"T1) . 

Besides the fermionic-like solutions found above, the Dirac-Kahler equation have also 
bosonic-like solutions that can be found by using an odd number of nodes in the con- 
struction of Dg and to impose 27r-periodic conditions on them. In this case, the 
eigenvalues are found to be ±\Jk{k + 1), k — 0, 1, • • - , with multiplicity 2{2k + 1) for 
k 7^ and 6 for k = 0. Since the method yields stable eigenvalues up to 16 digits of 
precision, i.e. they do not change as Ng and N v grows, the eigensolution in this case 
should be given in terms of trigonometric polynomials in 8 and (p. 



5 Final remarks 

The differentiation matrices introduced here to obtain the discrete form of an elliptic 
operator on a compact manifold have proved to be useful in exploring its spectral problem. 
This is not an unexpected feature. Examples of compact manifolds are spheres and tori, 
and on these, some parametrizations involve periodic or antiperiodic conditions on the 
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eigenfunctions which is assumed to be smooth enough in order to be expanded in a Fourier 
series or approximated by a trigonometric polynomial. Since the N x N differentiation 
matrix ([T]) yields the exact derivative of a trigonometric polynomial of degree (N — l)/2, 
a great number of nodes will produce a sufficiently approximate result which becomes the 
exact solution whenever this solution can be written in terms of trigonometric polynomials. 
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